EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 




CERN-PH-EP-20 12-147 
17 May 2012 



Experimental investigation of transverse spin asymmetries in /i-p SIDIS 

processes: Sivers asymmetries 



The COMPASS Collaboration 



Abstract 



The COMPASS Collaboration at CERN has measured the transverse spin azimuthal asymmetry of 
charged hadrons produced in semi-inclusive deep inelastic scattering using a 160 GeV beam and 
a transversely polarised NH3 target. The Sivers asymmetry of the proton has been extracted in the 
Bjorken x range 0.003 < x < 0.7. The new measurements have small statistical and systematic un- 
certainties of a few percent and confirm with considerably better accuracy the previous COMPASS 
measurement. The Sivers asymmetry is found to be compatible with zero for negative hadrons and 
positive for positive hadrons, a clear indication of a spin-orbit coupling of quarks in a transversely 
polarised proton. As compared to measurements at lower energy, a smaller Sivers asymmetry for 
positive hadrons is found in the region x > 0.03. The asymmetry is different from zero and positive 
also in the low x region, where sea-quarks dominate. The kinematic dependence of the asymmetry 
has also been investigated and results are given for various intervals of hadron and virtual photon 
fractional energy. In contrast to the case of the Collins asymmetry, the results on the Sivers asymme- 
try suggest a strong dependence on the four-momentum transfer to the nucleon, in agreement with 
the most recent calculations. 
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In the late 60's a simple and powerful description was proposed for the nucleon as a stream of partons 
each carrying a fraction x of the nucleon momentum in a frame where the nucleon momentum is infinitely 
large. From the dependence of the deep inelastic lepton-nucleon scattering (DIS) cross section on the 
energy and momentum transfered to the nucleon it was possible to identify charged partons with the 
earlier postulated quarks, and assess the existence of gluons as carriers of half of the proton momentum. 

Since the 90's it is well known that in order to fully specify the quark structure of the nucleon at twist- 
two level in quantum chromodynamics (QCD) three types of parton distribution functions (PDFs) are 
required: the momentum distributions q(x) (or ff(x)), the helicity distributions Aq(x) (or g\ (x)) and the 
transversity distributions Ajq(x) (or h\(x)), where x is the Bjorken variable. For a given quark flavour q, 
q(x) is the number density, Aq(x) is the difference between the number densities of quarks with helicity 
equal or opposite to that of the nucleon for a nucleon polarised longitudinally, i.e. along its direction 
of motion, and the transversity distribution Ajq{x) is the corresponding quantity for a transversely po- 
larised nucleon. If the quarks are assumed to be collinear with the parent nucleon, i.e. neglecting the 
intrinsic quark transverse momentum kj, or after integration over kj, the three distributions q(x), Aq(x) 
and Ajq{x) exhaust the information on the internal dynamics of the nucleon. On the other hand, from 
the measured azimuthal asymmetries of hadrons produced in unpolarised semi-inclusive deep inelastic 
scattering (SIDIS) and Drell-Yan (DY) processes a sizeable transverse momentum of quarks was derived. 
Taking into account a finite intrinsic transverse momentum kj, in total eight transverse momentum de- 
pendent (TMD) distribution functions are required to fully describe the nucleon at leading twist QJJ] . 
Presently, PDFs that describe non-perturbative properties of hadrons are not yet calculable in QCD from 
first principles, but they can already be computed in lattice QCD. In the SIDIS cross section they appear 
convoluted with fragmentation functions (FFs) iSS], so that they can be extracted from the data. 

A TMD PDF of particular interest is the Sivers function Aq q (or f^j), which arises from a correlation 
between the transverse momentum kj of an unpolarised quark in a transversely polarised nucleon and 
the nucleon polarisation vector |]4j|. In SIDIS this kj dependence gives rise to the "Sivers asymmetry" 
Ash, which is the amplitude of the sinO^ modulation in the distribution of the produced hadrons. Here 
the azimuthal angle is defined as = tp/, — (j> s with tpt, and <p s respectively the azimuthal angles of 
hadron transverse momentum and nucleon spin vector, in a reference system in which the z axis is the 
virtual photon direction and the xz plane is the lepton scattering plane. Neglecting the hadron transverse 
momentum with respect to the direction of the fragmenting quark, the Sivers asymmetry can be written 
as 



Z q e 2 q -q®D h q 



(1) 



where (£> indicates the convolutions over transverse momenta, e q is the quark charge and D h q describes 
the fragmentation of a quark q into a hadron h. 

In the very recent years, much attention has been devoted to the Sivers function, which was originally 
proposed to explain the large single-spin asymmetries observed in hadron-hadron scattering. The Sivers 
function is T-odd, namely it changes sign under naive time reversal, which is defined as usual time 
reversal but without interchange of initial and final state. For a long time the Sivers function and the 
corresponding asymmetry were believed to vanish [5] due to T-invariance arguments. However Brodsky 
et al. [a] showed by an explicit model calculation that final-state interactions in SIDIS arising from gluon 
exchange between the struck quark and the nucleon remnant (or initial state in DY) produce a non-zero 
asymmetry. One of the main theoretical achievements of the recent years was the discovery that the 
Wilson-line structure of parton distributions, which is necessary to enforce gauge invariance of QCD, 
provides the possibility for non-zero T-odd transverse momentum dependent (TMD) PDFs. According 
to factorisation the T-odd PDFs are not universal. The Sivers function can be different from zero but 
must have opposite sign in SIDIS and DY H. A lot of interest in the Sivers function arises also from 
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its relation with orbital motion of quarks inside a transversely polarised nucleon. In particular it was 
shown that orbital angular momentum must exist if the Sivers function doesn't vanish. Even though 
no exact relation between Sivers function and orbital angular momentum was derived yet, work is going 
on, also because the importance of assessin g th e role of the orbital angular momentum in the nucleon 
spin sum rule has grown in time (see e.g. [§-11]). 



Presently, the measurement of the Sivers asymmetry in SIDIS is the only direct way to assess the Sivers 
function. It became an important part of the experimental programs of the HERMES and COMPASS 



experiments, and it will be an important part of future SIDIS experiments at JLabl2 [ 12]. Furthermore, 



in the near future several experiments using the DY process will address the Sivers function, in particular 



its sign, in order to establish the prediction of restricted universality fll3Lll4ll . 



Using a 160 GeV longitudinally polarised /i + beam COMPASS measured SIDIS on a transversely po- 
larised deuteron ( 6 LiD) target in 2002, 2003 and 2004. In those data no sizeable Sivers asymmetry was 



observed within the accuracy of the measurements 1115141711 . a fact which is understood in terms of a 
cancellation between the contributions of u- and d-quarks. By scattering the e~ and e + beams at HERA 
off a transversely polarised proton target, HERMES measured in 2004 a non-zero Sivers asymmetry for 
positively charged hadrons [18]. A combined analysis of the COMPASS and HERMES data allowed for 
a first extraction of the Sivers function for u- and d-quarks Ill9l42lh . Still, as in the case of the Collins 
asymmetry, measurements on protons at higher beam energies were needed to disentangle possible higher 
twist effects. 

In 2007 COMPASS measured for the first time SIDIS on a transversely polarised proton (NH3) target. 
The results II22|] on the Sivers asymmetry for positive hadrons were found to be different from zero and 
turned out to be somewhat smaller than the final HERMES data H23TI . However the COMPASS results 
had larger statistical errors and a non-negligible overall scale uncertainty of ±0.01. A more precise 
measurement was thus mandatory and the entire 2010 data taking period was dedicated to this purpose. 

In this Letter, the results of the 2010 run are presented. They confirm with considerably smaller un- 
certainties the observation of the 2007 measurements. The higher statistics allow for first studies of the 
kinematic dependence of the asymmetry in a domain larger than the usual COMPASS DIS phase space. 

The COMPASS spectrometer is in operation in the SPS North Area of CERN since 2002. The principle of 
the measurement and the data analysis were already described in refs. 1 15- 17, 22, 24)]. The information 



on the 2010 run, the amount of data collected, the event reconstruction and selection, the statistics of 



the final samples, are given in a parallel paper on the Collins asymmetry B25I1 that was measured using 
the same data. In order to ensure a DIS regime, only events with photon virtuality Q 2 > 1 (GeV/c) 2 , 
fractional energy of the virtual photon 0.1 < y < 0.9, and mass of the hadronic final state system W > 5 
GeV/c are considered. A charged hadron is required to have at least 0. 1 GeV/c transverse momentum 
Pj with respect to the virtual-photon direction and a fraction of the available energy z > 0.2. This is 
refered to as "standard sample" in the following. 

The Collins and Sivers asymmetries are the amplitudes of 2 of the 8 azimuthal modulations, which are 
theoretically expected to be present in the SIDIS cross section for a transversely polarised target. They 
are extracted simultaneously from the same data as explained in ref. 12511 . The measured amplitude of 
the modulation in sin<I>s is £5 = /PrAsiv where / is the dilution factor of the NH3 material, and Pj 
the magnitude of the proton polarisation. In order to extract A$i V , the measured amplitudes £5 in each 
period are divided by / and Pj. The dilution factor of the ammonia target is calculated for semi-inclusive 
reactions [26] and is evaluated in each x bin; it increases with x from 0.14 to 0.17, and it is assumed 
constant in z and p\. The proton target polarisation (~ 0.8) was measured individually for each cell and 
each period. The results for A$i v from all periods of data taking are found to be statistically compatible 
and the final asymmetries are obtained by averaging the results from the full available statistics. Extensive 
studies were performed in order to assess the systematic uncertainties of the measured asymmetries, and 
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it was found that the largest contribution is due to residual acceptance variations within the data taking 
periods. In order to quantify these effects, various types of false asymmetries are calculated from the final 
data sample assuming wrong sign polarisation for the target cells. Moreover, the physical asymmetries 
are extracted splitting the events according to the detection of the scattered muon in the spectrometer (top 
vs bottom, left vs right). The differences between these physical asymmetries and the false asymmetries 
are used to quantify the overall systematic point-to-point uncertainties, which are evaluated to be 0.5 
times the statistical uncertainties. The only relevant systematic scale uncertainty, which arises from the 
measurement of the target polarisation, is evaluated to be 3% of the target polarisation. 

Figure Q] shows the Sivers asymmetries for positive and negative hadrons extracted from the 2010 proton 
data as a function of x, z and p\ , where the other two variables are integrated over. For negative hadrons 
the asymmetry is compatible with zero, while for positive hadrons it is definitely positive and stays 
positive down to x ~ 10~ 3 , in the region of the quark sea. There is good agreement with the published 
results from the COMPASS 2007 run Il22ll but with a considerable reduction of more than a factor of 
two in the statistical and in the point-to-point systematic uncertainties. Also, the asymmetry for positive 
hadrons is clearly smaller than the corresponding one measured by HERMES 12311 . This fact persists 
even when considering only events with x > 0.032, in the same x range as the HERMES experiment. 
The asymmetries in this restricted x range are shown as open points in fig. [2] 

The correlation between the Collins and the Sivers azimuthal modulations introduced by the non-uniform 
azimuthal acceptance of the apparatus as well as the correlations between the Sivers asymmetries mea- 
sured when binning the same data alternatively in x, z or p\ were already given in ref. 02511 . All correla- 
tion coefficients are found to be smaller than 0.2 and are relevant only in case of simultaneous fits of the 
various asymmetries. 

In order to further investigate the kinematic dependence of the Sivers asymmetry and to understand 
the reason of the difference with HERMES, the kinematic domain is enlarged to examine the events 
with smaller y values (in the interval 0.05 < y < 0.1), which correspond to smaller Q 2 and W values. 
Additionally, the standard data sample is divided into two parts, corresponding to0.1<y<0.2 and 
0.2 < y < 0.9. Since at small y there are no low-* data, only events with x > 0.032 are used. Figure [3] 
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Fig. 1: Sivers asymmetry as a function of x, z and p\ for positive (top) and negative (bottom) hadrons. 
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Fig. 2: Sivers asymmetry as a function of z and p\ for positive (top) and negative (bottom) hadrons. The 
open points (o , slightly shifted horizontally) are the values obtained in the range 0.032 < x < 0.70. The 
closed points (•) refer to the full x range and are the same as in fig-U 
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Fig. 3: Sivers asymmetry as a function of x, z and p\ for positive (top) and negative (bottom) hadrons 
for x > 0.032 in the y bins 0.05 < y < 0.1 (closed squares, ■), 0.1 < y < 0.2 (open triangles, V, slightly 
shifted horizontally) and 0.2 < y < 0.9 (open squares, □) . 



shows the Sivers asymmetries measured in these three bins of y as a function of x, z, and pj respectively. 
No particular trend is observed in the case of the asymmetries for negative hadrons (bottom plots), which 
stay compatible with zero as for the standard sample. A clear increase of the Sivers asymmetry for 
positive hadrons is visible for the low-y data. This strong effect can not be due to the slightly different 
mean values of x, since the Sivers asymmetry does not exhibit an x dependence for x > 0.032. On 
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Fig. 4: Left panel: mean value of y vs W. Middle panel: mean values of W vs x for the standard sample 
0.1 <y < 0.9 (closed circles, •) and for the samples 0.05 < y < 0.1 (closed squares, ■), 0.1 < y < 0.2 
(open triangles, V), and 0.2 < y < 0.9 (open squares, n). Right panel: mean values of Q 2 vs x for the 
standard sample 0.1 < y < 0.9 (closed circles, •) and for the samples 0.05 < y < 0.1 (closed squares, ■), 
0.1 < y < 0.2 (open triangles, V), and 0.2 < y < 0.9 (open squares, □). 
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Fig. 5: Comparison between the measured and calculated Sivers asymmetries for positive hadrons as a 
function of z for 0.1 < y < 0.9. The closed points (•) refer to the full x range and the open points (o) to 
the 0.032 < x < 0.70 range. The curves are from ref. [29]. 



the contrary, it could be associated with the smaller values of Q 2 and/or with the smaller values of the 
invariant mass of the hadronic system W. A similar dependence of the asymmetries on y was already 
noticed in the published results from the 2007 data. As can be seen from fig. |4] (left panel), there is a 
strong correlation between the y and W mean values: the mean values of W in the high x bins are about 
3 GeV/c 2 for the sample 0.05 < y < 0. 1 and larger than 5 GeV/c 2 for the standard sample 0.1 < y < 0.9 
(middle panel of fig. [4]). On the other hand, as can be seen in the right panel of fig. [4] bins at smaller 
y have smaller values of (Q 2 ). In particular, in each x bin the Q 2 mean value decreases by about a 
factor of 3 for the sample 0.05 < y < 0.1 with respect to the standard sample. Although the situation 
might be different in the target fragmentation region 027|], in the current fragmentation region the Sivers 
asymmetry is not expected to depend on y (or on W), while some Q 2 dependence should exist due to the 
Q 2 evolution of both the FFs and the TMD PDFs. 

Very recently first attempts to estimate the impact of the Q 2 evolution of the Sivers function l28ll led 
to encouraging results. In ref. [29] the Sivers asymmetry was evaluated for the HERMES kinematic 
region using the Sivers functions of ref. [30] and then evolved to the COMPASS kinematic region. The 
measured z dependence of the Sivers asymmetries for 0.1 < y < 0.9 is compared with the calculated 
one in fig. [5J for the entire x region and for x > 0.032. The linear trend of the data up to z — 0.75 is 
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Fig. 6: Sivers asymmetry as a function of x, z and pj. for positive (top) and negative (bottom) hadrons for 
0.032 < x < 0.70 in 3 different z bins: 0.1 < z < 0.2 (closed squares, ■), 0.2 < z < 0.35 (open triangles, 
V, slightly shifted horizontally when plottes vs x and p'j) and 0.35 < z < 1.0 (open squares, n). 



well reproduced, as well as the small increase of the slope for the high x sample. A very recent fit 113 ill 



of the HERMES asymmetries M23H and the COMPASS deuteron 111711 and proton [32] results given here 
was performed taking into account the Q 2 evolution in all x bins. It reproduces all the data well and 
provides strong support to the current TMD approach, which foresees a strong <2 2 -dependence of the 
Sivers function. 

We have also investigated the behaviour of the Sivers asymmetries at low z. Our standard hadron selection 
requires z > 0.2 to stay well separated from the target fragmentation region. In the range 0.1 < z < 0.2 
no effect on Astv is visible for negative hadrons, but one observes a clear decrease of the asymmetry for 
positive hadrons. In fig.[6]the data are plotted in 3 different z regions: 0.10 < z < 0.20, 0.20 < z < 0.35, 
and 0.35 < z < 1.00. While the shape of the asymmetry as a function of x stays the same, the size of 
the asymmetry shows a clear proportionality with z, in qualitative agreement with the expected linear 



behaviour (see, e.g. Q33|]). 



All the results given in this Letter are available on HEPDATA 13411 . The asymmetries for the standard 
sample as functions of x, z and pf have also been combined with the already published results from the 
2007 run Q and are also available on HEPDATA. 

In summary, COMPASS has obtained precise results on the Sivers asymmetry in SIDIS using a polarised 
proton target. A first investigation of its dependence on various kinematic variables shows significant 
dependences on z and y. By now, the Sivers asymmetry for positive hadrons is shown to be different 
from zero in a broad kinematic range and to exhibit strong kinematic dependences. After two decades 
of speculations, this is an important new insight into the partonic structure of the nucleon. In the light 
of the most recent theoretical advances refined combined analyses to evaluate the Sivers function and 
its dependence on the SIDIS variables are required in order to understand the role of the Sivers function 
in the various transverse spin phenomena observed in hadron-hadron collisions and in future Drell-Yan 
measurements. 
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